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A USER'S GUIDE TO THE SUDAN COMPUTER 
PROGRAM FOR DETERMINING THE VIBRATION 
MODES OF STRUCTURAL SYSTEMS 

By 

Raymond G. Kvaternik and Barbara J. Durllng 
SUMMARY 

The use of the SUDAN computer program for analyzing structural systems 
for their natural modes and frequencies of vibration Is described. SUDAN 
Is based on a substructures methodology and Is Intended for structures which 
can be represented as an equivalent system of beam, spring, and rigid-body 
substructures. A finite element stiffness technique In combination with a 
lumped raabs method which allows for non-diagonal masses Is used to model 
each of the substructures which comprise the system. User-written constraint 
equations are used to analytically Join the mass and stiffness matrices of 
the substructures to form the mass and stiffness matrices of the complete 
structure from which all the frequencies and modes of the system are 
determined. A feature of the SUDAN program Is the capability to treat the 
case In which both the mass and stiffness matrices of the coupled system may 
be singular simultaneously. A general description of the computer program 
Is given, the computer hardware and software specifications are Indicated, 
and the Input required by the program Is described. The usage of the program 
is illustrated by a simple example. The SUDAN program and related documentation 
can te .r ained from the Computer Software Management and Information Center 
(COSM",,) at the University of Georgia. 
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INTRODUCTION 


Two computational procedures for analyzing complex structural systems for 
their natural modsr and frequencies of vibration using a substructures method- 
ology are presented in reference 1. One of the procedures described is a 
direct method based on solving the eigenvalue problem associated with a finite- 
element representation of the complete structure. This computational procedure 
has been Implemented in a special-purpose computer program designated SUDAN 
(SUbstructurlng in Direct ANalysls). The SUDAN program is intended for struc- 
tures which can be represented as an equivalent system of beam, spring, and 
rigid-body substructures. A finite-element stiffness technique in combination 
with a lumped mass method which allows for non-diagonal masses is used to 
generate the mass and stiffness matrices for each of the substructures com- 
prising the system. The mass and stiffness matrices of the individual sub- 
structures are Joined analytically by means of user-written equations of 
constraint to form a finite-element model of the complete structure from which 
all the frequencies and modes of the system can be determined. A feature of 
the SUDAN program is the capability to treat the case in which both the mass 
and stiffness matrices of the coupled system may be singular simultaneously. 

The SUDAN computer program is somewhat limited in the maximum size of the 
problem which can be accommodated since it employs an in-core eigensolutlon 
technique. However, the program has been used quite successfully by the 
authors and their colleagues in several airframe vibration analyses carrying 
engineering responsibility. This experience has demonstrated a wide range 
of engineering applicability for the program. 


The purpose of this report is to present a user's guide for the SUDAN 
program. A general description of the program is given, the computer hard- 
ware and software specifications are indicated, and the inputs required by 
the program are described. A sample case consisting of an assembly of beams 
configured in the shape of an airplane and supported by springs is provided 
to illustrate the manner of preparing input data and to indicate the format 
of the output. The SUDAN program and related documentation can be obtained 
from the Computer Software Management and Information Center (COSMIC) at the 
University of Georgia. 

The present report relies on reference 1 for a description of the. analy- 
tical basis of the computational procedures which are. implemented in SUDAN 
and on reference 2 for a discussion of representing a structure as an assembly 
of beam, spring, and rigid-body substructures. For these reasons, references 
1 and 2 and the present report should he regarded as companion reports. 

GENERAL DESCRIPTION OF PROGRAM 

The SUDAN program is a system comprised of two programs — DSTIFF and BJD5 — 
which are executed sequentially with disk files providing the communication 
between them. 

Organization 

Based on a user-established idealization of a structure as an equivalent 
system of beam, spring, and rigid body substructures. Program DSTIFF employs 
the combined lumped-mass/finite element stiffness technique described in 
reference 2 to generate the stiffness and mass matrices for the structural 
block(s) comprising each of the constituent substructures. The uncoupled 
system stiffness and mass matrices, given by the composite matrices containing 
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the stiffness and mass matrices of the Individual substructures as submatrices 
on the principal diagonals, are then formed. Spring and/or mass terms coupling 
either the blocks of a substructure submatrix or t. . substructure submatrices, 
as well as any other additional spring or mass terms, are then added to the 
uncoupled system stiffness and mass matrices. User-written constraint equa- 
tions enforcing Inter-substructure displacement compatibility are then applied 
according to the method of Walton ar J Steevea (ref. 3)* leading to the gener- 
alized eigenvalue problem equation 

[K] {X} - LAMBDA [m] tx) (1) 

where [K] Is the coupled system stiffness matrix, [M] Is the coupled 
system mass matrix, LAMBDA Is the eigenvalue which corresponds to the square of 
the natural frequency, and {X> is the eigenvector which corresponds to the 
mode shape In terms of generalized coordinates arising from solution of the 
constraint equations. 

Program BJD5 solves the generalized eigenvalue problem given by equation 1. 
This program was written by the second author and is based on an algorithm 
devised by Walton and Durling (unpublished work, 1966) to reduce the free vib- 
ration equations of motion to symmetric standard eigenvalue form in the general 
case when both the stiffness and mass matrices may be singular. The analytical 
basis of the BJD5 algorithm is reviewed in reference 1. The solution provides 
all the eigenvalues and the eigenvectors of equation 1. Since the solution 
eigenvectors are generalized mode shapes, they are transformed back to the 
physical coordinates of the original uncoupled system for interpretation by 
the user. 

*See Appendix of reference 1 for description of an efficient computational 
procedure for implementing the method of Walton and Steeves. 
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continuity between reference* 1 and 2 and the present report Is facilitated 
by the use of the sam? notation whenever possible and through extensive use of 
comments In the source programs. 

Computer Hardware and Software Specifications 

The SUDAN programs — DSTIFF/BJD5 — and their subroutines are written In 
FORTRAN IV. They have been complied using the FORTRAN (FTN) compiler and 
executed on the Network Operating System (NOS 1.2) on Control Data Corporation 
(CDC) CYBER 173/175 digital computers at the Langley Research Center (LaRC) 
Computer Complex. These programs were developed on CDC 6400/6600 computers 
using the RUN compiler on a modified SCOPE 2.3 Operating SYSTEM (1COPS) at 
the LaRC Computer Complex. The only program change required to convert the 
programs from ICOPS to NOS was one end-of-flle test In Program DSTIFF. Library 
subroutines JACTV and MATINV are included In the source programs. The 
Intrinsic functions ABS and MOD and the basic external function SQRT 
utilized by the programs and subroutines are presumed to be available on other 
digital computer systems. Experience gained by the second author In converting 
programs received from other sites indicates that only a few changes would be 
required to make these programs operational on other digital computers. 

Compatibility features .- Some of the FORTRAN IV features employed or 
avoided in SUDAN in an effort to provide compatibility among various computers 
are as follows: 

1. All variables and arraya are Identified by a symbolic name consisting 
of from one to six letters or digits beginning with a letter. 

2. All variables and arrays are typed implicitly by the first character 

of the name: I, J, K, L, M, and N imply type INTEGER, and the other letters 
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Imply type REAL. No explicit type statements are used. 

3. Only returns of the type RETURN are used in the subroutines. 

When a RETURN statement Is encountered control Is transferred to the next 
statement following the CALL statement in the calling program. 

4. There are no ENTRY statements in the subroutines. 

5. Arrays have either one or two dimensions. 

6. Only one 2-branch IF statement is used in SUDAN, and it is the 
statement containing the end-of-file test in Program DSTIFF. 

7. Adjustable dimensions are formal arguments in all of the subroutines. 

8. The H specification is used in FORMAT statements containing a string 
of alphanumeric or special characters including blanks rather than delimiting 
the string with special symbols. 

9. Real constants formerly contained in arithmetic and relational 
expressions in a given routine have been replaced by real variables. These 
variables are not Included in a blank COMMON or a labeled COMMON statement. 

In each routine these variables are assigned values in DnTA statements and are 
not redefined during execution. Each constant in the DATA statements contains 
from two to fourteen digits, a decimal point, and an exponent comprised of the 
letter E followed by a signed or an unsigned integer constant. The three 
constants containing fourteen digits are in DATA statements in Subroutines 
JACTV and FREQ. 

10. Hollerith constants use the H specification. The letter H is 
preceded by the unsigned Integer 4 and is followed by four alphanumeric or 
special characters including blanks. In these programs, Hollerith constants 
are used only to define elements of the HMTX array. Since the HMTX array is 
included in a labeled COMMON statement, and since some computers require that a 


BLOCK DATA subprogram be used to Initialize arrays or variables in COMMON with 
DATA statements, DATA statements are not used to define elements of the HMTX 
array. 

Conversion considerations . - It is reconinended that the SUDAN programs and 
subroutines be converted to double precision If the computer to be used 
provides less than approximately fourteen digits In single precision computa- 
tion. All real variables and arrays except PROBID and HMTX in Program BJD5 
should be double precision. As a rule of thumb, when the difference of 
exponents (spread) between the largest eigenvalue and the smallest non-zero 
eigenvalue for a given eigenvalue problem approaches the number of digits of 
precision provided by the computer, the number of significant digits in the 
smallest eigenvalue may decrease to the point where even the first digit is not 
correct. 

Some specific modifications and points to be considered in converting the 
SUDAN programs and subroutines for use on other than CDC computers are as 
follows : 

1. Delete the PROGRAM statement in both Programs DSTIFF and BJD5 to make 
them main programs, and provide file assignments for the programs as required 
on a given computer. 

2. The precision of real variables and real arrays In the routines can 
be extended to double precision by Including the appropriate type statement 
for the given computer. The type statements may be DOUBLE PRECISION, IMPLICIT, 
or REAL with the extended length specification or combinations thereof. The 
precision of REAL constants can be increased by changing the letter E in the 
exponent of each constant to the letter D. REAL constants requiring 
increased precision occur only In DATA statements. It is noted that the 
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precision of two arrays — HMTX and PROBID— should not be extended slnas these 
arrays are used for Hollerith data. The arrays and the routines in which 
they are used are as follows: 

HMTX - DSTIFF, WKTXC, r . D5, PMTXC, FREQ 
P ROB ID - DSTIFF, BJD5 

3. If extended precision is provided in the routines, the appropriate 
number of digits should be added to the 14-digit constants which are in 
DATA statements in Subroutines JACTV and FREQ. 

4. A variable named ONEMIL is initialized to 10 b in DATA statements 
in both DSTIFF and BJD5, and it is not redefined during execution. This 
variable is used in determining the number of positive, non-zero eigenvalues 
since "zero" eigenvalues may not be precisely zero in the computer. Experience 
has indicated that 10 6 is an appropriate value for this test on CDC computers, 
but the value may have to be changed on other computers. 

5. The value range of a single precision (real) word on CDC computers 

-293 +322 

is 10 to 10 . The values of two variables, RHO and R1 , used in the 

programs are system dependent. The variable RHO In Subroutine JACTV specifies 

the accuracy requirement, and it is initialized to .IE-13 in the DATA 

statement. The variable R1 and its reciprocal R2 in Subroutine MATINV 

are used to keep the magnitude of the computational values within the range 

allowed on a given computer. The variable R1 Is initialized in the DATA 
100 

statement to 10 and is not redefined during execution. A programmer at 
a given site should be able to provide appropriate values for both RHO and 
Rl. 

6. If extended precision is provided in the routines, change each 
occurrence of the intrinsic function ABS to DABS and of the external 
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function SQKT to DSQRT. The functions and routines In which they occur 
are listed with the numbers enclosed In parentheses Indicating the number 
of occurrences of the function in the routine, as follows: 

ABS - JACTV(2) , FREQ(l) , MATINV(IO) 

SQRT - JACTV(l), BJD5(1) , FREQ(l) 

7. Program DSTIFF contains a two-branch IF statement which Includes 
an end-of-flle test (EOF function). This statement, which Is numbered 251, 
will have to be deleted on some systems and the appropriate end-of-f lie-test 
parameter Inserted In the READ statement, numbered 250, which immediately 
precedes it. 

8. Formatted input and output statements and the routines In which they 
occur are as follows: 

READ f ,11st - DSTIFF 

PRINT f - DSTIFF, BJD5 

PRIN: f ,11st - DSTIFF, BJD5, PMTXC, FREQ 

WRITE(6,f)llst - WMTXC 

where f represents the FORMAT number, list represents the Input/output 
(I/O) list to be transmitted, and 6 represents the output unit designator. 

In the PROGRAM statements for DSTIFF and BJD5, TAPES is equivalenced to INPUT 
and TAPE6 is equivalenced to OUTPUT. Although READ and PRINT statements are 
known to be valid FORTRAN formatted I/O statements on two non-CDC computers, 
they may have to be changed before the programs can be used on some computers. 

Problem size limitations .- The dimensions of the arrays In the source 
programs reflect a tailoring to the size of the sample problem. Larger 
problems can be accommodated by changing dimensions of the arrays and 
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Che values of Che varlaUes sec equal Co Che various dimensions and recom- 
piling Che programs. This 1» noc a formidable Cask since all of Che 
subroudnes utilize variable dimensions and che program source decks concaln 
derailed commentary. If Che value of Che variable MAXORD and dimensions 
assoclaced with lc are Increased In Program DSTIFF, Che user Is caudoned 
Co make appropriate changes In Program BJD5. B utilizing che maximum 
central memory (330K octal words) available Co a Job on Che CDC CYBER 173 
computers aC Che LaRC Computer Complex, the user can use approximately 
two hundred discrete coordinates Co define che uncoupled system. 

The primary reason for che problem size limitation Is that Che SUDAN 
programs use an in-core eigensoludon technique — designated Subroutine JACTV. 
This subroutine incorporates che threshold variation version of che serial 
Jacobi algorithm. Since this method is particularly appropriate for the 
special requirements of the SUDAN programs, the user Is cautioned against 
substituting another algorithm without first understanding the reasons for 
selection of JACTV as detailed in reference 1. If a replacement algorithm 
is used, It is recommended that comparison runs be made in both versions of 
the programs for a number of different problems to ensure that the Integrity 
of the programs has been maintained. 

General Input/Output Information 

Before attempting to prepare data for a calculation, it is advisable to 
read through this report and also to be familiar with the substructure 
modeling techniques and the procedure for formulating the equations of 
constraint which are described in reference 2. 

Only Program DSTIFF requires card input. Each of the programs produces 
both printed (formatted) and sequential binary (unformatted) output. The 
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binary output of Program DSTIFF ia tha binary Input for Program BJD5. Tha 
binary output of Program BJD5 may ba used, for example, as Input to a user- 
selected program for ganaratlng and annotating mode-shape plots. In sub- 
scquant sections of this paper, tha Input and the output of tha programs 
are described, and the card Input for and the printed output from executions 
of a sample case are listed. 

The card Input data for Program DSTIFF must be prepared using the pre- 
scribed formats and must be assembled In tha proper sequence. A consistent 
system of units must be used in preparing the applicable data. Integers 
must be Input with correct values since many of the Integers, or variables 
computed from them, are used to control the number of values read Into 
arrays. The data described In the tables are designated either as required 
(REQ) or as lonal (COND). As implied by the designator, all required 

data must be prepared. If conditional data are Included, the appropriate 
lnteger(s) must be input as Indicated in the tabular description of these 
data. The number of cards required for each read statement is tabulated 
or can be computed from the expression which is tabulated in the column 
v Jer "No. of Cards." Note that the quotient for each of these Integer 
expressions must be incremented by one if the dividend is not an exact 
multiple of the divisor. 

The identification card should contain alphanumeric Information which 
will be meaningful to the user for later identification of the case. Integers 
must always be right Justified in the specified fields. More latitude Is 
allowed for the single-precision, floating-point numbers: a decimal point 

in the field overrides the decimal position specified in the format, and an 
exponent is not required if a properly positioned decimal point is Included. 
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Whan an exponent is Included, It oust be right justified In th<> field with 
the letter E and/or the sign of the exponent preceding it. 

As noted previously, the binary output of Program OSTIFF is the binary 
input for Program BJD5. The data are contained on two files — TAPE 3 and 

TAPE4. No file manipulation Is required between executions when these 
programs are executed sequentially on the same job. When the programs are 
executed on separate Job submittals, the user must save both TAPE 3 and TAPE 4 
after the execution of Program DSTIFF on the first Job, and then must retrieve 
these files on the second Job prior to the execution of Program BJD5. The 
file names, TAPE3 and APE4, are declared, respectively, on the DSTIFF pro- 
gram statement as the third and fourth files and on the BJDS program state- 
ment as the second and third files. In order to avoid errors on file 
assignments, the user is cautioned that Program BJD5 uses TAPE4 as a scratch 
flic for temporary data storage after it reads the coupled mass and stiffness 
matrices from it. Program BJD5 performs no write operations on TAPE3. 

The binary output from Program BJD5 is contained on TAPE9. If this 
file is needed for a later job submittal, it must be saved by the user prior 
to the termination of the Job on which it is created. The file name TAPE9 
is declared as the fifth file on the BJDS program statement. 

DESCRIPTION OF SUDAN INPUT AND OUTPUT 

The SUDAN computer program system ij comprised of two programs — 

DSTIFF and BJDS — and their subroutines. Usually these programs are executed 
sequentially with disk communication between them. Program DSTIFF calls 
subroutines JACTV, WMTXC, and ZEROM: Program BJD5 calls subroutines FREQ, 
JACTV, MATINV, and PMTXC. The subroutines are Included in the SUDAN source 
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program. A sample deck, structure for execution ot the SUDAN programs on the 


CVBER 173/175 computers at the iaRC Computer Complex Is 

This section Is divided Into four subsections: 

1. Input data for Program DSTIFF. 

2. Output for Program DSTIFF. 

3. Input for Program BJD5. 

4. Output for Program BJD5. 

Input Data for Program DSTIFF 

The data deck for Program DSTIFF may consist of as 
types of data sets Identified and assembled as follows: 

1. Initial data 

2. Block data 

3. Additional mass data 

4. Additional stiffness data 

5. Constraint data 
Initial data, block data, and constraint data must always be prepared; the 
additional mass data and/or the additional stiffness data are Included when 
appropriate for the case. Definition of a beam substructure requires at 
least one block of data and may require up to four blocks (bending In two 
principal directions, torsion, and extension), depending on the freedoms 
assigned to the beam component. It Is noted, for clarity, that the data 
deck for Program DSTIFF Is contained in one (1) record. 


shown in figure 1. 



many as five different 
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Initial data.- The initial data consist of two cards. 


No. of 
Cards 


Option FORMAT 


Description of Initial Data 


1 REQ 20A4 PROBID • one card containing alphanumeric 

Identification for the case (maximum 
of 80 characters) 


1 REQ 614 NBLKS , NORDER, NCEQ8, NMASAD, NSPRAD, IPRINT 

NBLKS - number of blocks in the uncoupled 
system stiffness and mass matrices 

NORDER - order of the uncoupled system 

N'lEQS »- number of constraint equations 

NMASAD » number of mass terms to be added 
on or above the diagonal of the 
uncoupled system mass matrix. The 
program will provide symmetry. 

NSPRAD - number of spring terms to be added 
o n or above the diagonal of the 
uncoupled system stiffness matrix. 
The program will provide symmetry. 

IPRINT ■ control integer for printing output 

IPRINT ■ 1, standard output- will be 
printed for the case 
IPRINT - 2, the assembled uncoupled 

stiffness and mass matrices 
as well as the standard 
output will be printed. The 
order of these two square 
matrices is NORDER. The 
additional printing may be 
helpful, especially to a new 
user, in determining that the 
conditional additional mass 
and/or stiffness terms have 
been input correctly. 


The two initial data cards are followed by NBLKS sets of block data. 

Block data .- There are NBLKS sets of block data. Each set provides the 
input data for generation of one block of the uncoupled system stiffness and 
mass matrices. The first card in each set contains five integers. 
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No. of 
Cards 


Option FORMAT 


Input Variables with Description for 
Card 1 of a Block Data Set 


1 REQ 514 KK, LOC, MK, MM, ISYM 

KK - order of the block which ia equal to the 
number of degrees of freedom of the 
block 

LOC - control Integer for block type 

LOC « 1, beam bending 

LOC - 2, beam torsion 

LOC ■ 3, beam extension 

LOC « 4, spring or rigid body component 

MK - control integer for spring components 

MK - 1, read stiffness matrix from cards 
MK ■ 0, no read, block will be null 

MM - control Integer for rigid body components 

MM ■ 1, read mass matrix from cards 
MM ■ 0, no read, block will be null 

ISYM ■ control Integer indicating whether 
symmetry is used in formulating the 
problem 

ISYM - 1, the program will multiply the 
block stiffness and mass 
matrices by 2.0. This option 
is used only when a component 
is out of the plane of sym- 
metry and the complete problem 
is being solved as a pair of 
symmetric and antisymmetric 
problems 

ISYM ■ 0, no multiplication - either not 
applicable to component or the 
complete problem is being 
solved without taking advan- 
tage of symmetry. 

The user completes preparation of a given set of block data from the table 
wh^jh describes the additional input variables required for the block type 
specified by the value of LOC. The data descriptions for both beam torsion 
(LOC - 2) and beam extension (LOC « 3) are included in the same table since 
they are similar. 


15 


No. of Option FORMAT Additional Input Variable* with Description for 

Cards Beam Heading, LOC ■ 1, Block Data 

1 REQ 314 X, N, KASE 

The free-free stiffness matrix is generated in 
the partitioned form 

I "I 
A | B 

B* J C 

where A is KxK, B is KxN, B 1 is NxK, and 
C is NxN. For the free-free case, K ■ N » 
(KK)/2 

K • order of matrix A 
N ■ order of matrix C 

KASE ■ control integer for mass matrix 

KASE ■ 1, diagonal mass only 
KASE * 2, diagonal mass and rotary inertia 
KASE “ 3, non-diagonal mass and rotary inertia. 

The off-diagonal terms for this block 
must be input under the NMASAD option. 


K/5 

REQ 

5E14.8 

X array of K values, the running local coordinate 
positions (dimensional) 

(K-l)/5 

REQ 

5E14.8 

S array of K-l values, the average stiffness, El, 
between a pair of sequential coordinates 

K/5 

REQ 

5E14.8 

DM array of K values, lumped masses located at 
each of the K coordinates (zero values are 
admissable) 

M/5 

COND 

5E14.8 

RI array of N values, lumped rotary inertia for 


each of the N rotational coordinates (zero 
values are admissable) . These values must be 
input if KASE is not equal to 1. 


If the bean bending data do not complete the NBLKS-th set of block data, the 
user must turn back to the beginning of the block data section to start the 
next set of block data for the case. If, however, the set number for beam 
bending data is equal to NBLKS, data preparation is continued with the next 
type of data which may be conditional additional mass data, conditional 
additional stiffness data, or constraint data, in that order. 
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No. of 
Cards 


1 


Additional Input Variables with Description 
Option FORMAT for Beam Torsion, LOG - 2, and 

for fleam Extension. LOC - 3. Block Data 

REQ 214 K, KASE 

K - order of free-free stiffness matrix and 
beam torsional inertia or extenslonal 
mass matrix 

KASE - identification integer for mass 
matrix 

KASE - 4, diagonal torsional inertia or 
extenslonal mass 

KASE - 5, non-diagonal torsional inertia 
or extenslonal mass. The off- 
diagonal terms for this block 
must be input under the NMASAD 

option. 


K/5 

REQ 

5E14.8 

X array of K values, the running local 
coordinate position (dimensional) 

,i’-l)/5 

REQ 

5E14.8 

S array of K-l 

values 




S array - GJ 
S array * AE 

for beam torsion (LOC-2) 
for beam extension (LOC- 

K/5 

REQ 

5E14.8 

RI array of K 

values (zero vrlues are 


admissible) 

RI array • lumped torsional Inertias for 
beam torsion (LOC-2) 

RI array « lumped extenslonal masses for 
beam extension (LOC-3) 


If the beam torsion or beam extension data do not complete the NBLKS-th set 
of block data, the user must turn back to the beginning of the block data 
section to start the next set of block data for the case. If, however, the 
set number for beam torsion or beam extension data is equal to NBLKS, data 
preparation is continued with the next type of data which may be conditional 
additional mass data, conditional additional stiffness data, or constraint 
data, in that order. 
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No. of Option FORMAT Additional Input Variables with Description 

Cards for Spring or Rigid-Body Components, 

1.0C - 4, Block Data 

(KK> ”- / 5 COND 5E14.8 A array of KK x KK values, the stiffness 

matrix for a spring component. It Is read In 
continuously by rows If the value of MK for 
the block is not zero. Each card has five 
values except the last card which may have 
blank fields to the right of the last value. 

The program provides a null mass matrix for a 
spring component. 

(KK) 4 '/5 COND 5E14.8 B array of KK x KK values, the mass matrix for 

a rigid-body component. It Is read In 
continuously by rows if the value of MM for 
the block is not zero. Each card has five 
values except the last card which may have 
blank fields to the right of the last value. 

The program provides a null stiffness matrix 
for a rigid-body component. — NOTE — The 
case in which the rigid body has a null mass 
matrix Is admlssable. This option is selected 
by setting both MK and MM to zero. 

If the spring or the rigid-body component data do not complete the NBLKS-th 
set of block data, the user must turn back to the beginning of the block data 
section to start the next set of block data for the case. If, however, the set 
number for the spring or the rigid body component data Is equal to NBI.KS, data 
preparation Is continued with the next type of data which may be conditional 
additional mass data, conditional additional stiffness data, or constraint data, 
in that order. 
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A dditional mass data .- If the centers of gravity of the lumped masses 
for a beam substructure are not located on the elastic axis there will be mass 
terms coupling the blocks In the mass matrix corresponding to the beam 
substructure. Additional mass terms (diagonal, of '-diagonal , and coupling) 
will also arise If th u mass matrix of a rigid body component not treated as a 
substructure Is combined with the mass matrix of the beam component to which 
It Is attached. The use of consistent mass matrices for bending, torsion, or 
extension will lead to non-diagonal mass matrices for those blocks. The user 
Is referred to reference 2 for discussion of these aspects. The control 
Integer NMASAD determines whether additional mass data are to be read. All 
off-diagonal mass data for beam components are Input as additional mass data. 
The row number and the column number provided with each of the NMASAD terms 
are the row and the column of the element in the uncoupled system mass matrix 
to which the additional term Is to be added. The program will provide symmetry, 
so data are Input only for values which appear on or above the diagonal. 

CAUTION: For a symmetric or an antl-symmetrlc formulation cf a problem, 

the masses to be added to the uncoupled system mass matrix must be multiplied 
by 2.0 by the user prior to input If they correspond to components out of the 
plane of symmetry. 


No. of 
Cards 


Option 


FORMAT 


Input Variables with Description for 
Additional Mass Data 


NMASAD/4 COND 4(213, E14. 8) Row number, column number, and value for each 

of the NMASAD additional mass terms to be 
added. Blank fields may not be present In 
any except the last card, and then may be 
only those fields to the right of the 
NMASAD-th group of three numbers. Note 
caution in tue previous paragraph. 

Data preparation is continued with conditional additional stiffness data If 


NSPRAD is not zero and with constraint data If NSPRAD equals zero. 


1 ? 



Additional stiffness data. - Spring terms which couple the substructure 


submatrices will arise If the stiffness of springs having neither end tied to 
ground are combined with the stiffness matrices of the components to which they 
are attached. For springs which have one end tied to ground, additional spring 
terms (diagonal, off diagonal, and coupling) will also arise If the other end 
of the spring is attached to a beam or rigid body component at a point other 
than a station established as a degree of freedom In the problem formulation. 

The user Is referred to reference 2 for discussion of these aspects. The 
control integer NSPRAD ditermines whether additional stiffness data are to be 
read. The row number and the column number provided with each of the NSPRAD 
terms are the row and the column of the element In the uncoupled system stiff- 
ness matrix to which the additional term Is to be added. The program will 
provide symmetry, so data are input only for values which appear on or above 
the diagonal. 

CAUTION: For a symmetric or an anti-symmetric formulation of a problem, 

the spring terras to be added to the uncoupled system stiffness matrix must be 
multiplied oy 2.0 by the use r prior to input if they correspond to components 
out of the plane of symmetry. 

No. of Option FORMAT Input Variables with Description for 

Cards Additional Stiffness Data 

NSPRAD/4 COND 4(213, E14. 8) Row number, column number, and value for 

each of the NSPRAD additional spring terms 
to be added. Blank fields may not be 
present in any except the last card, and 
then may be only those fields to the right 
of the NSPRAD-th group of three numbers. 
Note caution in the previous paragraph. 

Data preparation is continued with constraint data. 
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t oust r.i im il.it a . - Since all elements of the constraint matrix are set to 


zero before constraint data are re’d, onlv non-zero elements are Input. The 
number of rows In the constraint matrix Is equal to the number of constraint 
equations, NCEQS. The number of columns in the constraint matrix Is equal to 
the order of the uncoupled system, NORDER. In the description which follows, 
NZECM is used only to indicate the total number of non-^ero elements of the 
constraint matrix and Is neither an Input parameter nor a program variable. 


Option Input v. tri.ii>;, th Description 

Cards Constraint Data 


NZECM/4 REQ 4(213, E14. 8) Row number, column number, and coefficient 

for each non-zero value In the constraint 
matrix. Blank fields may not be present 
In any except the last card, and then may 
only be those fields to the right of the 
last (NZECM-th) group of three numbers. 

The maximum value for a row number is 
NCEQS and for a column number is NORDER. 


Output for Program DSTIFF 

P rinted outpu t. - The printed output for Program DSTIFF consists of a 
block by block summary of the input data arranged in tabular form and the 
eigenvalues of the constraint eigenvalue problem which is associated with the 
method of Walton and Steeves (ref. 3) for imposing the equations of constraint. 
An option is provided for printing the uncoupled system mass and stiffness 
matrices. This additional output may be helpful in identifying input errors 
in the placement of the additional mass and stiffness elements in the uncoupled 
system mass and stiffness matrices. The printed sample output presented in a 
subsequent section does not include this optional output. 

Binary output . - The binary output of Program DSTIFF is not described 
here since it is fully described as the binary input for Program BJI)5. 
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Input for Program BJD5 


Card Input .- There le no card Input for Program BJD5. 

Binary Input .- The binary Input for Program BJD5 la the binary output 
from Program PSTIFF, and It la contained on two files, TAPE3 and TAPE*. 

Thla Input lnciudea problem Identification, numbers which specify the order 
of matrices contained on the files, the elements of the connectivity matrix 
(BETA) which enforces displacement compatibility at the Interfaces of the 
substructures, and the elements of the coupled system mass and stiffness 
matrices. The BETA matrix Is used at the end of Program BJD5 to transform 
the modal matrix back to the original physical coordinates of the uncoupled 
system to facilitate Interpretation by the user. 


File Record 

Name Number 


Description of the Binary Input for Program BJD5 


TAPE3 1 


2 through 
(NRBETA + 1) 


PROBID, NRBETA, NCBETA 

PROBID • the card Image (80 characters) of the 

problem Identification card provided as 
input for Program DSTIFF 

NRBETA • number of rows In the BETA matrix 

NCBETA - number of columns in the BETA matrix and 
order of both the coupled system mass 
matrix and the coupled system stiffness 
matrix 

NRBETA records each of which contains NCBETA elements 
for a row of the BETA matrix 


TAPE 4 


1 This record contains all the elements of the coupled 
system mass matrix of order NCBETA. 

2 This record contains all the elements of the coupled 
system stiffness matrix of order NCBETA. 
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Output tor Program BJD5 


Printe d outpu t. - The printed output for Program BJD5 conalsta of the 
eigenvalues of the coupled system maaa matrix, the eigenvalues and frequencies 
of the coupled system, and the eigenvectors of the coupled system transformed 
back to the original physical coordinates of the uncoupled system. 

Binary output .- The binary output of Program BJD5 la contained on one 
file, TAPE9, and Includes the number of rows and columns In the modal matrix, 
problem identification, frequencies, and vectors of the modal matrix which 
have been transformed back to the original physical coordinates of the 
uncoupled system. These data may be used, for example. In a program selected 
by the user to generate and annotate mode-shape plots. If this output Is 
required for a later Job submittal, It must be saved ji her on a magnetic 
tape or on a permanent file prior to termination of the Job on which It la 
generated. 

CAUTION: In order to provide compatibility among various computers, 

the PROBID array contains the 60 characters of the problem Identification 
card In 20 words each of which hs* 4 left-justified characters. When the 

a 

PROBID array Is to be used to annotate the plot(s) generated by a user- 
supplied plotting program, the user is reminded that the 80 characters are 
contiguous only If the computer word has four 6-blt bytes (32 bits). If, 
for example, the computer word has 36 bits or 60 bits, several suggestions 
to assist the user are presented, as follows: 

1. Plot four characters at a time from each of the 20 words In the array 

2. Reformat the 80 characters In memory from the 20-word array to a 
compatible array for the given computer. For example, assuming PROBID is 
dimensioned 20 and CASE Is dimensioned 8, the transfer could be accomplished 
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on CDC computers by Including the following two statements In the program: 
ENCODE (80, 5000, CASE) PROBID 
5000 FORMAT (20A4) 

Since the word size on CDC computers Is 60 bits, the eight words In the CASE 
array will contain the 80 contiguous characters after execution of the ENCODE 

statement. 

3. Modify Program BJD5 so that the 80 characters are reformatted Into 
contiguous locations In an array as described In the second suggestion, and 
then write the new array Instead of PKOBID on TAPL9. 

4. Perform a binary (unformatted) read of only NRB and K from the first 
record on TAPE9, and provide an Input problem Identification data card to be 
read by a formatted read statement which will place the 80 characters in 
contiguous locations in an array on the given computer. Ou CDC computers, 
the format would be 8A10, and the array would be dimensioned 8. 


File Record 

Name Number 


% 

Description of the Binary Output from Program BJD5 


TAPE9 1 NRB, K, PROBID 

NRB - number of rows in the modal matrix after it has 
been transformed back to the original physical 
coordinates of the uncoupled system 

K - number of frequencies and also number of columns in 
the modal matrix 

PROBID * the card image of the problem identification 
card provided as input for Program DSTIFF. 

It is noted that the 80 characters are con- 
tained in 20 words each of which has 4 left- 
justified characters. See caution which pre- 
cedes this table. 

2 K calculated frequencies (Hz) of the coupled system 

arranged in ascending order. 


3 K records each of which contains the NRB elements of a 

through vector (column) of the modal matrix after it has been trans- 

(R + 2) formed back to the original physical coordinates of the 

uncoupled system. The k-th vector is associated with the 
k-th frequency. 
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SAMPLE CASE 10 ILLUSTRATE USAGE 01 PROGRAM 


Description of Conf lguratlon 

The SUDAN program waa applied to the vlbratl< \ analysis of a model con- 
sisting of an assembly of beams configured In lu shape of an airplane in 
reference 1. This model is the basis of The b<> pie -:ase selected to illustrate 
the usage of the program. The geometric properties of the mouel are given in 
figure 2. Since attention will be directed to a symmetric vibration analysis 
as In reference 1, advantage can be taken of symmetry and only one hr If of 
the model need be explicitly considered. Thus, for analysis purposes, the 
mathematical model can be considered to ^ mist of three substructures: 
fuselage, wing, and tall. For symmetric motions, the component motions of 
Interest are vertical bending of the fuseluge beam and vertical bending and 
torsion of the wing and tail beams. The distributed mass of the fuselage, 
wing, and tail beams is lumped at discrete points along the elastic axes of 
the respective members such that the sectional mass, mass unbalance about the 
elastic axis, and the torsional inertia are preserved. Although the rotary 
inertia in bending associated with each ol th< oe lumped masses is not zero, it 
is assumed unimportant an. taken to be zero. Since the mass axis and elastic 
axis of each beam comprising the model are coincident, the static unbalance of 
the wing and tall masses about their respe-otivr elastic axes is zero. However, 
to illustrate the use of the NMASA1) option in the program, the static unbal- 
ances of the wing masses about the elastic axis are assumed to be non-zero. 

This will inertially couple wing bending with wing torsion. A sketch of the 
resulting lumped-mass model of the beam assembly is given in figure 3. The 
fuselage and wing support springs shown in figure 3 were not part of the 
configuration analyzed In reference 1 but are introduced here to Illustrate 
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Che use of Che NSPRAD option In Che program. The fuselage aupporc springs 
consisC of a Cransladonal spring having a spring race and a rocadonal 

spring having spring race k , which are accached Co Che elasdc axis aC Che 
sevench fuselage sCadon. The wing sutporc springs are cransladonal springs 
having a spring race k^ and are accached Co Che leading edge of Che wing 
(2.54 cm from elasdc axis) aC Che local axis sCadon corresponding Co Che 
fourch wing station. Since Chese springs are accached Co Che wing off Che 
Che elasdc axis, additional coupling of wing bending and wing torsion will 
oc"ur Chrough stiffness terms associated wlch Che wing supporC springs. 

Each fuselage station has two degrees of freedom: vertical translation 

and rotation. Each wing and Call station has three degrees of freedom: 
vertical translation, rotation, and torsion. These discrete degrees of 
freedom are Identified and ordered In Table I. A total of 65 discrete degrees 
of freedom are associated with the uncoupled system, 26 of which correspond to 
coordinates with no masj as a consequence of assuming the rotary Inertias to 
be zero. The model structural properties, a discretized for the analysis In 

accordance with the lumped-mass model shown In figure 3, are summarized in 

a 

Table II. The procedure for determining the numerical values of the stiffness 
coupling terms associated with the support springs and for identifying the 
position of the mass static unbalance terms and spring terms In the substruc- 
ture mass and stiffness matrices is described In reference Both the wing 

mass static unbalance terms and the stiffnesses of the support springs shown 
in Table II were deliberately chosen to be small so as to represent only small 
perturbations from the configuration treated in reference 1, thereby 

*It should be noted that the structural properties of the wing and tall 
beams are based on cuts perpendicular to the elastic axes of those components. 


tacllltat lng a comparison with the results shown therein if desired. A 


summary of the constraint equations required to analytically join the 
substructures is given in Table 111. 

Listings of the Input and printed output for lire sample case are given in 
tire next two sections, respectively. The standard output option (IPRINT " 1) 
was selected for tire sample case and thus the optional output (the uncoupled 
system masa and stiffness matrices'! is not included in the sample output 
list lng. 

Ihe sample case considered has no unrestrained (i.e., rigid-oody) degrees 
of freedom, and thus none of the eigenvalues LAMBDA are zero. If a structure 
is unrestruined, the eigenvalues corresponding to the rigid-body modes should 
be zero. In practice, however, the eigenvalues which are calculated for the 
rigid-body modes may not be exactly zero. Because some of these "zero" eigen- 
values may be negative, Program BJD5 takes the square root of the absolute 
value of LAMBDA to obtain the nat*ir u i frequency in radians per second. 



27 







i n 

E 





(A 

i- 






E 

0 ) 

H- 





C 

* 4 -» 

o 





41 








OS 

«A 

X 





c 






t/> 


C 

i- 




(A 

L. 

4 ) 

4 -' 

E 



4 

a 

E 

<* 

4 



E 

C/l 

4 ) 

E 

0 ) 





H 


A 



rH 

rH 

4 ) 

-♦-* 


E 

E 

4 

C* 


C 

II 


as 

C 

c 

O 


O) 

0 ) 

0 ) 

o 

o 

i_ 

c* 

4 

13 

XI 

•— 


4 ) 

L. 

H 

a> 



4 -* 

N 


0 ) 

H 


•• 

1 

01 

Ui 

c 

.» 

■a 

T 3 

C 

c 

3 

•• 

ft 

T> 

T 3 

o 

0 

U. 

* 


< 



o 



O) 

a> 

c 

o> 

C -X 

C -* 

O -X 

C -¥ 

— o 

— u 

•• <J 

— o 

*o o 

•o o 

Ui o 

T 3 O 

C rH 

C -H 

L. rH 

C -H 

0 ) «Q 

41 _Q 

O XI 

41 X 3 

CD 

JO 

•— 

CO 


C 

O -X 

— o 
U) o 
L. *H 

o n 


o 

o 








o 

9 




O JN 

A- <0 


•4 




AJ 


^ * 

-4 f4 

<o cy 

9 

<y 

r* 

<r 

*4 

o 


0 * 

• • 

AJ 

9 

o 

9 

-o 

•4 

o 


O r4 

A- A- 

« A> 

O 

• 

AJ 

o 

• 

o 

r 

4 <y 

m o 

O O 

# 

A- 

•4 

-r 


o 

4 

• • 

•4 mm 

• • 

• 

AJ 

• 

• 

•# 

• 

a 










O 










o 










a 










a 



















o 

r 









9 

4 









A* 

o 


(f> K> 

•o -c 


m* 

1 A 



AJ 

=> 

«■ ® 

-4 -4 

9 cy 


<y 


o 

•4 

o 

%/s 

^ AJ 

• • 

AJ AJ 


o 

cy 


*4 

o 


O —4 

A* A* 

A- A- 

o 

• 

AJ 

o 

• 

o 

m 

m ® 

IA fA 

O O 

m 

A- 

*4 

m 

«o 

o 

O 

• • 

•4 *4 

• • 

• 

AJ 

• 

• 

<# 

• 


UJ 




O »4 

A • • 

• #4 | 





> # * 

# ^ 

-4 ® >r 




<p 

•# <# 

9 ^ *\ 




9 






m 

o ^ 

AJ 9 


•4 A 


m* 

m # 


AJ 

<y « 

AJ 

•4 O 



n 

o «r 

A 

— o 


# 

o 

• o 

o 

• o 



AJ 

A- O 

AJ 

o 



• 

AJ • 

• 

^ • 

o 






O 9 






9 -C 

9 





o -n 

»A 





O AJ 

AJ 





Aj O Aj 

O 




<© 

O O • 

9 




cy 

OOKN 

9 O 




IA 

O • A 

a # a 


*4 


•4 

• 1 *4 

4 4 • 

•# 

cy m 

■# 

H O 



Aj 

o ^ 

AJ 

•4 O 

«A A- O 

9 9 0 

a> 

• o 

m 

• o 

* 9 

IA O 

-4 

A* O 

-4 

A O 



• 

AJ • 

• 

<# • 

9 m o 

AJ IA U> 





Aj IA JA 




N l 
(N» ISJ O H . 


O H Ki ^ I 
(\j K t I" i 


«C 9 «* 
9 9 
N N •*( 
A* A* 9 

o o o 


9 

c 


9 

AJ 


o 

o 

o 

o 


9 ® 
o 9 

• >0 
a. o 


<£> 

9 


r4 O 
*4 O 

• o 
*r> o 

9 • 


«v ® 
9 a- 


<4 iTl l/> 

OCW> 
O O IT 

o o o 


9 «# 

Aj N 

O O 

0 0 9 
• 9 O O 



»/» 

I/* »T* i 

9 

ST 


ui O l 

x 


•o o 

-4 9 

o o 

-4 O 


ia i/n ./> tf> 


9 • 

o a* r* i 


O O O 
9 AJ 9 

MON 


099 
• 9 o 


> ia ja #a o O O O -« 


o o 
9 A- 


9 9 
AJ o 
^4 0 0 H 


O O *4 

4 r-« 

o o • 

o tr* 

• 9 


O m k 
N H O » 


oir/>( 

44Hi 
H • • 

ONM 


o *a 

9 cy 
A- w 

AJ -4 -4 

o o 
o o 
o o 
o © o 


o -• 


O • H H H H 


i 9 o» o 
* 0-0 m 
ia aj o < 
IOiTNO 


o <y <y « a 
<roo 9 9 a«- a- 
O o • • O Aj 

• /» K N O H 

o t N N • • 


*> O 

99-40 
m a* 

•4 AJ 

OhhOO 
® -4 *4 O O AJ O 
O O • • O O <4 

• a> k> ✓> o o 

o • 9 9 • • 


O' ® 
o 9 
• 9 
A» O 


■a 

<y 


• o 

• o 


NOO 
9 9 AJ 


r n o 

AJ IA Al 


o 9 9 9 

^ O AJAJ JJ O 
O I A • fl W 
• AJ A- O O 

O • AJ • • 


• O 
•A O 


or a. 
cy 9 
O o 
o o 

0-400 
>r *4 o o 

O A • O O 

• AJ tf\ o O 

O • 9 • • 


9 «4 9 tf> 
* *> AJ 


*4 IA rf> O 


O AJ 

9 ® 

® 9 
® 9 

o 9 o 

O O • rf> 
o o o 9 
O O *A • 

• • m f 

*4 *> A. 9 

•# •# 9 

A* -4 A» © 
AJ IA IA 


o «r* 

• • 
*4 | 


A* # 
AJ <A 


O 
o « 

• o 


AJ o 

•4 O 


28 




Sample Printed Output for SUDAN 
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Last mode (corresponds to largest eigenvalue] 
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CONCLUDING REMARKS 


A user's guide to the* SUDAN computer . rogram for analysing Htructural 
systems for their natural modes and frequencies of vibration has been given. 
Sll, JN 1h baaed on a substructures methodology and la intended for atructurea 
which can be repreaented aa an equivalent ayatero of beam, spring, and rigid- 
body substructures!. A finite-element atlffneaa technique in combination with 
a lumped mass method which allows for non-diagonal masses is used to model 
each of the substructures which comprise the system. User-vritten constraint 
equations are used to analytically join the muss and stiffness matrices of 
the substructures to form the mass and stiffness matrices of the complete 
structure from which all the frequencies and modes of the system are 
determined. A general description of the program was given, the computer 
hardware and software specifications were indicated, and the input required 
by the program was described. The usage of the program was illustrated by a 
simple example. 
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TABLE I 


- IDENTIFICATION OF DECREES 01 FREEDOM FOR SYMMETRIC VIBRATION 
ANALYSIS OF AIRPLANE BEAM ASSEMBLY 


FUSE LACE 

q l " q 13 
q 14 ‘ q 26 

U'INC 

q 27 ' q 33 
q 34 ~ q 40 
q 41 " q 47 

TAIL 

q 48 " q 53 
q 54 " q 59 
q 60 " q b5 



Displacements 

Slopes 


Displacements 

Slopes 

Twists 


Displacements 

Slopes 

Twists 


} 

} 

} 

} 


Vertical bonding 


Vertical bending 
Torsion 


Vertical bending 
Torsion 





TABLE II.- DISCRETIZATION EMPLOYED FOR A ! RIMJINE BEAM ASSEMBLY 


Station 

Local 

Coordinate 

Foaltlon 

Maas* 

Torsional 

Inertia 

El 

CJ 

Mass 

Static 

Unbalance 


m 

k| 

k» 

2 

-m 

N-m 2 

N- 

2 

•m 

kg-m 

Fuselage 









1 

0.0000 

0.03648 

Not 

137.15 

Not 

0.0 

2 

0. 1016 

0.07296 

applicable 

137.15 

app 1 1 

cab U» 

0.0 

3 

0.2032 

0.07296 



137.15 



0.0 

4 

0. 3048 

0.07296 



137.1') 



0.0 

5 

0.4064 

0.06567 



137.15 



0.0 

fefc 

0.4877 

0.52653 



137.15 



0.0 

7 

0. 6096 

0.08026 



137.15 



0.0 

8 

0.7112 

0.07296 



137.15 



0.0 

9 

0.8128 

0.07296 



137.15 



0.0 

10 

0.9144 

0.07296 



137.15 



0.0 

11 

1.0160 

0.06567 



137.15 



0.0 

12 

1.0973 

0.07296 


f 

137.15 

\ 

f 

0.0 

13 

1.2192 

0.04378 


r 

■ 


r 

0.0 

Wing 









1 

0.0000 

0. 06*87 

0. 1395 

x 10 -4 

27.091 

45 

113 

+0.04448 x 10" 3 

2 

0. 1016 

0. 12975 

0.2790 


27.091 

45 

.113 

+0.08896 

3 

0.2032 

0.12975 

0.2790 


27.091 

45 

.113 

+0.13345 

4 

0. 3048 

0.12975 

0.2790 


27.091 

45 

.113 

+0.17793 

5 

0.4064 

0.12975 

0.2790 


27.091 

45 

.113 

+0.22241 

6 

0.5080 

0. 12975 

0.2790 


27.091 

45 

.113 

9-0.26689 

7 

0.6096 

0.06487 

0.1395 


■■ 

— 

-0.13345 

Tall 









1 

0.0000 

0.03243 

0.06979 x 10~ 4 

27.09) 

45 

.113 

0.0 

2 

0.0508 

0.06487 

0.13958 

27.091 

45 

.113 

0.0 

3 

0. 1016 

0.06487 

0.13958 

27.091 

45 

.113 

0.0 

4 

0. 1524 

0.06487 

0.13958 

27.091 

45 

.113 

0.0 

5 

0.2032 

0.06487 

0.13958 

27.091 

45 

.113 

0.0 

6 

0.2540 

0.03243 

0.06979 

_ 



0.0 


‘*Rotary Inertia of lulled masses neglected 

'includes mass of stem and coll of shaker used In reference 1 


kj - 350. N/m, k, ■ .0565 N-m/radlan, k^ - 87.6 N/m 
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TABLE III. 


‘ CONSTRAINT EQUATIONS FOR AIRPLANE 


BEAM ASSEMBLY 


lm q 27 " ^6 ‘ 0 

2 ’ q 34 Cos 0^+ q^ 8ln e . o 

3. -<I 34 sin 6j+ q 41 cos - 0 

4> q 48 - q i 2 * 0 

5 - ^54 cot* J 2 + q 6Q sine,- 0 

6 * “ q 54 sln 0 2 + q 6£) cos0 2+ q^ . o 


0 1 " 8W **P of wing elastic axis - 30° 
°2 - sweep of tail elastic axis - in 0 


Data for DSTirr 


knd-of-rarord card 


Directives for DSTirr 


/ 7/0/9 End-of-rtcord card 


COMflD n. Sava TAPK9 If aodae ara naadad for a later job aubatttal 


100 . 


LDsn i Pfttsrr a-m; i Nr . krk- fatal > 


mi( l«C0NPlLK.R0i]VD) 


*>Dinr(r.i) 


U0TIOPL-UD5/U0-49064 7M) 


HKW l MP (OPL . COMP I LE .LOO) 


LOO. 


ldskt t rRKsETA-w; t nr . im- fatal ) 


m(i»coMriLl,Koiris) 


».KT (OPL-DST I rr /UN-4906 A 7*0 



Figure 1 # - Sample deck structure for execution of SUDAN program on CYBER 173/175 
computers at the LaRC Computer Complex. 
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rh 


0.318 cm thick' 


All material steel 


i 30° /£ / 


5.08 cm 


V 

7 *> 


0.052 cm square • 


t 


48.77 cm 


53 


121.02 cm 



hl8ure 1 '~ Geometric properties of airplane bean assembly. 
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I 


*-Po , nt of attachment of springs kj and Ic, 


Figure 3 


(b) Planform view 


- Lumped-mass model of airplane beam assembly supported by springs. 
(Note: Offsets of masses and spring attachment points from wing 
elastic axis are exaggerated for illustrative purposes) 
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